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We analyze the Stern-Gerlach experiment in phase space with the help of the matrix Wigner 
function, which includes the spin degree of freedom. Such analysis allows for an intuitive visualization 
of the quantum dynamics of the apparatus. We include the interaction with the environment, 
as described by the Caldeira-Leggett model. The diagonal terms of the matrix provide us with 
information about the two components of the state, that arise from interaction with the magnetic 
field gradient. In particular, from the marginals of these components, we obtain an analytical 
formula for the position and momentum probability distributions in presence of decoherence, that 
show a diffusive behavior for large values of the decoherence parameter. These features limit the 
dynamics of the present model. We also observe the decay of the non-diagonal terms with time, and 
use this fact to quantify the amount of decoherence, from the norm of those terms in phase space. 

From here, we can define a decoherence time scale, which differs from previous results that make 
use of the same model. 


I. INTRODUCTION 

The Stern-Gerlach (SG) experiment is a cornerstone in 
quantum mechanics. It showed, for the first time, direct 
evidence for the discretization of the spin states of the 
electron, by analyzing the motion of Silver atoms through 
a magnetic field gradient [1]. Most textbooks make con¬ 
tinuous use of the SG, as a simple way to illustrate the 
quantum measurement process, since the electron spin 
only involves a two-dimensional Hilbert space. 

A consistent description of the SG experiment needs, 
obviously, to be quantum, even though one can make an 
introduction based on a semiclassical description, using 
a spin-dependent force that gives rise to a “trajectory” 
that depends on the initial spin state. The full quantum 
treatment reveals a richer dynamics, as it leads to entan¬ 
glement between the spin and spatial degrees of freedom 
[2-6]. 

In this paper, we perform a phase space analysis of the 
SG device, including the interaction with the environ¬ 
ment. This interaction will be described by the Caldeira- 
Leggett model [7]. In this respect, the starting point is 
similar to the analysis in [8]. However, our description 
is based on the use of the Wigner function (WF) [9]. 
Wigner functions have proven to be a powerful tool in 
physics, and can be used as an alternative formulation 
of quantum phenomena, including their dynamics. The 
particular features of the phase space description make it 
particularly advantageous in some situations, for instance 
recognizing the quantum features of states, or dealing 
with decoherence scenarios. In the WF, interference ef¬ 
fects manifest in a clear way [10-12]. 

In order to include the spin degree of freedom, one 
needs to extend the usual definition of the WF. A com¬ 
mon prescription in the literature is the use of a matrix 
valued WF [13], where the spin indices give rise to the 
matrix elements. Such description has some advantages 
when dealing with a particle subject to a spin-dependent 
force, since some effects like the spin precession, or mo¬ 


tion that depends on the spin component, are better vi¬ 
sualized with respect to a fixed spin basis. Examples 
of this description are a previous analysis of the Stern- 
Gerlach experiment without including decoherence ef¬ 
fects [ 14], the study of entangled vibronic quantum states 
of a trapped atom [13], or the reconstruction of the fully 
entangled quantum state for the cyclotron and spin de¬ 
grees of freedom of an electron in a Penning trap [15]. 

As we will show, the phase space description provides 
a clear visualization of the SG phenomenology. First, the 
diagonal terms show the motion of the two components 
of the quantum state (corresponding to spin up or down 
along the gradient direction). By considering an initial 
Gaussian state with arbitrary spin direction, we can ob¬ 
tain the marginals from the Wigner function, which de¬ 
scribe the appropriate probability distribution function 
(PDF) of position or momentum [16]. Each of these PDF 
have a Gaussian shape with a center and width which are 
modified by the interaction with the environment, and 
give valuable information about the state evolution. On 
the other hand, the out of diagonal terms can be used 
to describe the effect of decoherence, which manifest into 
a damping of the norm associated to these terms. We 
use this norm as a figure of merit to quantify the amount 
of decoherence experienced by the system, as a function 
of the parameter 7 that quantifies the strength of the 
coupling with the environment. As a result, we obtain a 
decoherence time which scales as 7 -1 / 5 , in contrast with 
previous results that claimed a 7 1 / 3 scaling [8]. We ar¬ 
gue that our result is more realistic, as it implies that 
a larger decoherence parameter manifests into a shorter 
decoherence time scale. 

The rest of this paper is organized as follows. In Sect. 
II, we solve the equations that govern the evolution of the 
matrix WF for the SG device, when the master equation 
based on the Caldeira-Leggett model is introduced to de¬ 
scribe the environment. By evaluating the marginals of 
the diagonal elements, we obtain the position and mo¬ 
mentum PDFs, and we analyze some limiting situations. 
Sect. Ill is devoted to the analysis of our results in a 
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standard setup of the SG experiment. In particular, we 
study the damping of the off-diagonal terms, and make 
use of this to define a decoherence time scale. Finally, we 
discuss the validity of the model to describe decoherence 
phenomena for the above setup. Our conclusions are pre¬ 
sented in Sect. IV, while some cumbersome expressions 
have been relegated to the Appendix. 

II. DYNAMICS OF MATRIX WIGNER 
FUNCTIONS INCLUDING INTERACTION 
WITH THE ENVIRONMENT 

The behavior of a spin 1/2 neutral particle under the 
action of a magnetic field gradient, including the influ¬ 
ence of decoherence effects, will be studied in this section. 
Particles entering the SG apparatus will move initially 
along the tube, defined as the x axis. The geometry of 
the magnetic field can be described by a dependence of 
the form 


{\z, a) = \z) (g> |a)} , where \z) Iz € ffi. is the eigenbasis of 
the position operator, and {|a)} is a fixed basis in spin 
space. We find it convenient to choose the eigenstates 
of S z (IS* = +h/2) = |+>, I S z = -h/2) = I-)). Finally, 
7 is the damping rate of the system in the environment. 
The coefficient D is defined as D = 2m r ykBT, with ks 
the Boltzmann’s constant, and T the temperature of the 
environment. 

As discussed in the Introduction, the analysis of the 
dynamics of this model will be presented on phase space, 
with the help of Wigner matrices 

1 ps s s 

W a p(z,p, t) = — J ds e~ lE ^(z+-,a\p{t)\z --,/?), 

^ (4) 

where W a p(z,p) are the spin matrix elements of the 
Wigner function in the above-mentioned S z base. The 
matrix WF has, among others, the following properties: 

1. One has 


B(x, y , z) = rjy] + ( B 0 - rjz)k, ( 1 ) 


which contains a uniform part B 0 , and a gradient of mag¬ 
nitude g on the plane orthogonal to the x axis. Notice 
that both position dependent terms in the latter equation 
are necessary in order to satisfy V • H = 0 and V x B = 0. 
However, it can be shown that the effect of the mag¬ 
netic held contribution along the y direction causes fast 
oscillations due to Larmor precession, which can be av¬ 
eraged out. Following [17], we neglect this contribution 
(see also [ 6 ]). In this way, in absence of decoherence ef¬ 
fects, the problem can be effectively factorized as the free 
propagation along the x and y axis, and the nontrivial 
motion corresponding to the z coordinate, which can be 
described by the Hamiltonian 


where p is the canonical conjugate momentum for z, and 
in and g s are the mass and gyromagnetic ratio of the 
particle, respectively. In Eq. (2), a z is the third Pauli 
matrix, and p,g the Bohr’s magneton. 

In order to account for decoherence effects, we assume 
that they are described by the Caldeira-Leggett master 
equation [18], which accounts for those effects in the sys¬ 
tem via collisions with a thermal bath of particles. In the 
position and spin representation, with the Hamiltonian 
( 2 ), the master equation can be written as [ 8 ] 
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In the latter equation, p a p{z, z ', t) = (z, a\ p(t ) \z', (3) are 
the matrix elements of the density operator p(t) repre¬ 
senting the particle state, at a given time t, on the basis 


Wp a (z,p,t) = W*p{z,p,t), (5) 

which implies that the matrix WF is Hermitian. 

2. The normalization condition becomes 

/ OO /»00 

/ W aa (z,p,t)dzdp = 1. ( 6 ) 

-oo J —oo 


3. The marginal distributions of (4) are related to ma¬ 
trix elements of the density operator. In particular, 
for the diagonal components we have 


/ OO 

W aa (z,p, t)dp = (z,a\p(t)\z,a) = f (±) (z,t), (7) 

-OO 


with a = dr, where f^(z,t) represents the posi¬ 
tion PDF for the particle. In a similar way, the 
marginal over the position variable 


W aa (z,p, t)dz = {p,a\p(t)\p,a) = g {±] [p,t) ( 8 ) 


represents the momentum PDF, with \p, a) = | p) <g) 
l a )' {|p} /p G R} being the eigenstates of the mo¬ 
mentum operator. 


Using the definition Eq. (4) and the dynamics of the 
density matrix Eq. (3), one can obtain the corresponding 
differential equations for the Wigner function, given as 
follows 


dW^\z,p,t) _ p dW { d k \z,p,t) +£) 9 2 H^ ±} ( 2 :,p,t) 


dt to dz 

OlyWpHz.V.t)) 

+7 - Si - 
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dp 2 


( 9 ) 
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dW od (z,p,t) _ p dW od (z,p,t) D d 2 W od (z,p,t) 
dt m dz dp 2 

, d(pW od (z,p,t)) 

+7 — oi — 

± 2i\{B 0 + i]z)Wqj(z,p, t) 
h 

where W d ± \z,p,t) stands for the diagonal elements of 
the Wigner function, with the upper sign corresponding 

to W++(z,p,t), and the lower sign to W ( z,p,t ). We 

also defined W od (z,p,t) = W H _( z,p,t ), which implies 

that W-+(z,p,t) = W* d (z,p,t), according to property 1 
above. 


A. General solution of the differential equations 

To solve the system of equations (9,10) a Fourier trans¬ 
form is performed over both z and p [8]. Once the equa¬ 
tions are solved, the Wigner functions are retrieved by 
performing the inverse Fourier transform. 

Assuming that the incident particle is described by a 
polarized Gaussian beam with spin | n) = a |+) + b |—) 
(with |o | 2 + | 6 | 2 = 1 ), the initial density operator p( 0 ) 
can be written as 

P(0) = | i>) (V’l ® \n) (n\ , (11) 

with the wave function that represents the initial state 
| - 0 ) in position space defined as 

' Kz) = (12) 

From the initial state Eq. (11) one can obtain the 
matrix Wigner function at t = 0, which can be written 
as 

W(z,p,0)=(\f b TjWi(z,p), (13) 


where 


0 ft . 2 0-2 

Wi(z,p) = - -- (14) 

7TO 

represents the Wigner function for a spinless Gaussian 
state, and a is the Gaussian width of the particle’s spatial 
PDF. 

Solving the differential equation (9) by the method 
commented above, with the help of the initial condition 
(14), one finds the general solution for the diagonal ele¬ 
ments. After some algebra, they can be written as follows 


W^iz^t) 


7I V<2(t) 


(15) 


and we have defined the new variable r = yt. The func¬ 
tions F(z,p, t) and G(r) are defined in the Appendix. In 
these functions we introduced the notations 


2c 


?yA (r + e T — 1) 

7 2 TO 


(16) 


Pc = 


p\ (1 — e T ) 

7 


(17) 


The role played by z c and p c will be discussed in the next 
Section. 

Using the same procedure for the differential equation 
( 10 ), the solution for the off-diagonal elements is also 
found. The resulting expression is lengthy so that, in 
order to express it in a more compact way, we introduced 
the functions C'j(r) (i = 1, 2, 3,4,5, 6 ), that can be found 
in the Appendix. The off-diagonal elements can finally 
be written as 


W od (z,p,t) 


-+c x - 
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where we have omitted, for simplicity, the dependence of C'j(r) on the variable r. The matrix Wigner function takes 
the following form 


W(z,p,t) 


( M 2 W^ + \z,p,t) ab* W od (z,p,t) \ 
\a*bW; d (z,p,t) \b \ 2 W d ~\z,p, t) J - 


(19) 


As we discuss below, the diagonal terms in W(z,p,t) describe the behavior of the particles in phase space, and the 
off-diagonal terms represent the coherence of the state. 


B. Marginals of the Wigner function: position and 
momentum PDFs 


In fact, it is referred to as a quasi-probability distribu¬ 
tion , and may even take negative values. This had to be 


The Wigner function (for a spinless particle) can not be 
associated with a probability distribution in phase space: 
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expected from first principles, given the incompatibility 
of the position and momentum observables in quantum 
mechanics. One can, however, obtain the PDF corre¬ 
sponding to the particle position by integrating over the 
momentum variable, and vice-versa, as described in the 
previous Sect. Eq. (7) can be integrated, with the result 

UT»c ) 2 

f {± \z,t)= 6 ——, (20) 

V7rcr z 


where 

2 _ 2D (2r + 4e~ r - e~ 2r - 3) H 2 (1 - e~ T f 2 

Gz q 3 m 2 + 7 2 m 2 a 2 ° 

( 21 ) 

is the squared width of the position distribution. Inte¬ 
gration in Eq. ( 8 ) leads to the momentum PDF: 


( p-Fpc) 2 

9 (±) {p,t ) = e —^ -, (22) 

V P 


with 


4 = 


2D (l — e~ 2r ) h 2 e~ 2T 


(23) 


giving the squared width of the momentum distribution. 

The above results for the marginals, Eqs. (20,22) 
clearly show that the diagonal components of the Wigner 
matrix correspond to Gaussian distributions in phase 
space ( z,p ) which center (±z c ,±p c ) and width depend 
both on time, and on the rest of parameters of the prob¬ 
lem, including the decoherence constants 7 and D. 

Let us notice the following properties of these quanti¬ 
ties: 


1. By differentiating Eqs. (16) and (17) one readily 
obtains 


dz c _ p c 
dt to ’ 


(24) 


dp c 

dt 


??A - 7 p c . 


(25) 


which can be easily identified as the classical equa¬ 
tions of motion for a particle subject to a constant 
force, plus a friction term. These equations allow 
us to describe the motion of the center of the two 
Gaussians using a semiclassical framework (espe¬ 
cially if we neglect the interaction with the envi¬ 
ronment, as done in most textbooks). 


2. Let us consider the limit yf <C 1. Performing a 
Taylor expansion gives 


. h H ,2Da 2 

p c ~rj\t, <jp — T (——- 1 ) 7 *. (27) 

00 7 n z 

If we further neglect the last term in Eq. (27), the 
above results can be easily interpreted as the action 
of a constant force on the particle, and agree with 
the ones expected for the experiment in a decoher¬ 
ence free environment [17]. 


3. In the opposite limit (i.e., when yt 1) we can 
approximate 


r?A + 

z c ~ - 1, o z 

7 m 



(28) 


Pc - —, Vp - 0 ™ = — = \J4mk B T (29) 

7 7 

which leads to a limiting value of the momentum 
center (and width as well). This is a well known ef¬ 
fect in classical mechanics, that appears under the 
action of a friction force, and will play an impor¬ 
tant role in our analysis when large values of 7 are 
involved. We also obtain o z oc y/t, i.e. the charac¬ 
teristic behavior of a diffusive regime. Fig. 1 fea¬ 
tures the evolution in time of the center of the PDF 
(both in position and momentum), using the values 
of the parameters as defined in the next Sect. For 
low values of 7 , the center z c of the position PDF 
first grows quadratically, and then it does linearly. 
The growth time scale is dictated by 7 _1 , so that 
for very large values of this parameter z c remains 
close to zero. On the other hand, the center p c of 
the momentum PDF grows linearly for moderate 
values of 7 , but approaches a constant value as 7 is 
increased, the asymptotic limit being inversely pro¬ 
portional to the decoherence parameter. In Fig. 2 
we have plotted the width of the position and mo¬ 
mentum distributions as a function of time. The 
position width o z grows quadratically for small val¬ 
ues of t (as compared to 7 -1 ), whereas the regime 
o z oc y/t corresponds to late times. The transition 
can only be seen for 7 = 1 on this Figure, given 
the 7 -1 scaling. As for the plots representing o p , 
one can check that the ratio 2D £> is of the order 
~ 10 13 (notice that this ratio is independent of 7 ). 
Eq. (27) then predicts a fast increase of o p even 
at early times, as it is clearly observed from these 
plots. This magnitude will then reach the asymp¬ 
totic value o also on the same 7 _1 time scale, 
which can only be appreciated, in that Figure, for 
the 7 = 10 10 case. 


III. APPLICATION TO THE STERN-GERLACH 
EXPERIMENT 


- 2 c 


r]Xt 2 
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The study of a realistic SG experiment setup with de- 
(26) coherence effects will be the core of this section. First, 
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Figure 1. (Color online) Plots of the center of the position 
(up) and momentum (down) PDFs, as given by Eqs. (16) 
(17), respectively, for three values of the parameter 7 . 


parameters for the setup will be introduced. Then, the 
evolution of the system will be pictured in phase space 
using the results from the previous section. To conclude, 
the characteristic decoherence time of the system will be 
studied through the damping of the off-diagonal elements 
of the Wigner function as the system evolves. 


A. Experiment parameters 

We assume an incident beam of Silver atoms (to = 
1.8 x 10 -25 kg, g s ~ 2) starting in the \S X = fi/ 2 ) = 
|+) + ^75 |—) spin state, with an average speed v = 500 
m/s, and a beam width a = 10 ^ 5 m. The SG apparatus 
parameters are based on the realistic ones used in a pre¬ 
vious work [19]. In this setup, the applied magnetic field 
is B 0 = 5 T, and the gradient 77 = 1000 T/m. The longi¬ 
tude of the tube is l = 0.2 m, which implies a flight time 
of around 0.4 ms for the above Silver atoms speed. This 
is, therefore, the characteristic time scale for the system 
dynamics. The operational temperature of the tube is 
around T — 300 K, i.e. the laboratory temperature. 




Figure 2. (Color online) Plots of the width of the position 
(up) and momentum (down) PDFs, for three values of the 
parameter 7 . 


B. Phase space representation 

In this section we will illustrate the behavior of the 
SG experiment in a decoherent environment, by showing 
some plots of the Wigner function using the parameters 
defined in Sect. Ill A. 


1. Diagonal elements 

In order to represent the phase space distribution of 
the particle, the trace of the matrix Wigner function is 
plotted for 7 = 1 , 10 3 , 10 10 s _1 in Figures 3, 4 and 5, 
respectively. The trace allows us to show the total quasi¬ 
probability distribution, thus putting on the same plot 
both spin components. Fig. 3 (7 = 1 s _1 ) shows how 
the incoming state splits into two separating components 
within the experiment timescale (< 1 ms): At t = 100 fj,s 
both terms start to split, and at t = 200 fj,s they are vis¬ 
ibly separated. We observe the distortion of the original 
shape of the Wigner function, caused by the different evo¬ 
lution in 2 and p, that appears even when the interaction 
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Figure 3. (Color online) Contour plots of the trace of the matrix Wigner function with 7 = 1 s 1 . The left panel corresponds 
to the initial (t = 0 ) state, while the middle panel shows the situation at t = 100 fis, and the right panel is for t = 200 p,s. 
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Figure 4. (Color online) Same as Fig. 3, for a value 7 = 10 3 s 1 


with the environment is not included (see [14]). Note 
also that friction with the environment quickly broadens 
the spatial width, from the initial value to the millime¬ 
ter scale in just t = 200 ps. The same behavior can be 
observed for the momentum width. 

As can be seen in Fig. 4 (7 = 10 3 s _1 ), a higher value 
of the friction force quickly limits the momentum of the 
peak centers, and causes the growth of the distribution 
widths. The influence of the SG apparatus is still visible 
in the shape of the beam, but due to the speed limit and 
the continuous growth of the width, the two components 
of the beam do not separate. An even larger value of the 
friction force (Fig. 5) causes the distribution peak centers 
to remain at the origin, and the width of the momentum 
distribution quickly grows, thus completely masking all 
the effects of the apparatus. 


2. Off-diagonal elements 

To study the off-diagonal Wigner function, two 3D- 
plots of the real part of W 0 d{z,p,t ), corresponding to 
7 = 1 , 10 10 s _1 , are drawn in order to see how the value 
of the damping constant affects the decoherence rate. At 
t = 0 (Fig. 6 ) this amounts to representing the initial 
Gaussian distribution Eq. (14). 

As time goes on, Figs. 7 and 8 show how the real part 


of W 0 d(z,p, t) evolves, from a Gaussian distribution, to 
an oscillatory function. We also observe that these os¬ 
cillations are damped due to the interaction with the 
environment. By comparing both figures, we immedi¬ 
ately see that a larger value of 7 increases the oscillation 
frequency and reduces the amplitude of the oscillations, 
thus leading to a faster decoherence. 


C. Decoherence time 

In order to quantify the loss of coherence in the sys¬ 
tem we choose, as a figure of merit, the norm A(t) of 
W a d{z,p,t) at a given t, defined as 

/ oo roo 

/ \W od (z,p,t)\ dz dp. (30) 

-00 «/ —OO 

This quantity provides the total volume, in phase space, 
occupied by the coherent term of the Wigner function. 
Taking the module, instead of the function itself, avoids 
for cancellations due to the oscillatory nature of this 
term. In Fig. 9 we have plotted A(t) for different values 
of 7 . As expected, decoherence effects manifest in a de¬ 
crease of A(f) with time. In fact, the approach to zero 
appears at earlier times as one considers larger values of 
7 , clearly indicating a faster decoherence process. 
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Figure 6 . (Color online) 3D-plot of the real part of W 0 d(z, p , t) 
at t = 0 . 
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Figure 8 . (Color online) Same as Fig. 7, for 7 = 10 10 s 1 



Re(W od ) 


p (g m/s/mol) 


z(pm) 


Figure 7. (Color online) 3D-plot of the real part of W 0 d(z, p, t ) 
at t = 0.005 ps (up) and t = O.lps (down), for 7 = 1 s -1 . 


In view of this result, we can introduce a decoherence 
time td, as the time it takes for A (t) to reduce its initial 
value by a factor e. From the above data, we can obtain 
td for a given value of 7 . These data are collected in Fig. 
10 . 



Figure 9. (Color online) Plot of A(f), as defined by Eq. (30), 
for different values 7 = 1 , 10 2 , 10 4 , 10 6 , 10 8 s _1 , as a function 
of time. 





















Figure 10. (Color online) Red points show the decoherence 
time td for the same values of 7 used in Fig. (9). The bar 
indicates the imprecision that originates from the data. The 
solid curve is our fit td = a 7 (see the text for explanation). 

Performing a numerical fit to the curve td{ 7 ) we find 
that the decoherence time can be approximated by the 
formula td = a 7 b , with a = 0.497 ± 0.002 and b = 
—0.198 ± 0.001, pointing to a behavior as td oc 7 -1 / 5 . 
This result is at variance with the one discussed in [ 8 ], 

( 3 h 2 m 2 2 \ 

4 £>™ 2 A 2 J was 

claimed, which translates into a dependence of the form 
td oc 7 1 / 3 , i.e. the larger the decoherence parameter 
7 , the later the system experiences decoherence. In our 
opinion, such result is unrealistic, as one expects the op¬ 
posite behavior: Decoherence should take place faster as 
7 is increased, in accordance to our results. 

D. Discussion about the validity of the model 

In the SG experiment, Silver atoms travel through the 
high vacuum beam pipe. We assume that they collide 
with the residual air molecules, resulting into Brownian 
motion and decoherence. In the Caldeira-Leggett master 
equation, these effects are included via the damping rate 
7 . The Langevin equation relates 7 to the viscosity ll as 
follows [ 20 ] 


with // the viscosity of the medium for Silver atoms. This 
quantity can be explicitly calculated through the follow¬ 
ing expression [ 21 ] 

/r = 67 t R\j! = 2.994 tt RnM mvm Am, (32) 

where R = 144 pm is the radius of the Silver atom, n 
is the particle number density of the air; Mm, Vm and 
A m are the molecular mass, mean speed and mean free 


path of air molecules, respectively. The Stokes’ law term 
67 tR accounts for the contribution of Silver atoms mass 
to the viscosity [22], and the term // = 0.499 uM m vm^m 
accounts for the air contribution. Assuming typical val¬ 
ues for the pressure of the beam pipe P ~ 10 -6 mbar 
[23], one finds [21] n ~ 10 10 cm -3 , Am — 10 4 cm, and 
vm — 500 m/s. To obtain a rough estimate, we assume 
that the molecules in the low pressure conditions in the 
tube are mostly composed by IV 2 , with a mass Mm — 28 
uma = 4.65 10 " 26 kg. 

Using the above data, one obtains the estimate /i ~ 
10 ~ 14 Kg/s, implying that the order of magnitude of 7 
is ~ 10 10 s _1 . Such value, however, posses a problem: 
As we have seen (Fig. 4), for damping rates larger than 
10 3 s _1 , the two components of the initial state do not 
separate due to the high friction with the environment. 
This is obviously not what happens in the real experi¬ 
ment, leading us to serious doubts of the application of 
the Caldeira-Leggett model to the SG experiment. 

A possible explanation of this disagreement might be 
that the obtained value of 7 is wrong, because the cal¬ 
culation of the viscosity discussed above is not correct in 
high vacuum. As [21] indicates, when Am is larger than 
the size of the container, the gas is in a molecular state, 
and cannot be characterized by a viscosity anymore. Fur¬ 
thermore, as explained in [23], due to the high vacuum 
in the beam pipe, the mean free path of Silver atoms is 
a multiple of the beam pipe length. Under these con¬ 
ditions, one should question the validity of the master 
equation Eq. (3), which was used to derive the corre¬ 
sponding differential equations that describe the Wigner 
function dynamics. 


IV. CONCLUSIONS 

In this paper, we have studied the SG experiment 
with environmental induced decoherence described by 
the Caldeira-Leggett model. Our description is done on 
the phase space, making use of Wigner functions, with 
the additional spin degree of freedom. Our goal was to 
describe the kinematics of the atoms traversing the mag¬ 
netic field gradient, and interacting with a thermal bath 
of particles, leading to Brownian motion and decoher¬ 
ence. We solved the differential equations for the Wigner 
function corresponding to the model, starting from an 
initial separable state, with a Gaussian shape in space, 
and an arbitrary spin state. The diagonal terms on the 
S z basis have a simple interpretation, in terms of the 
two separating states of the apparatus. By calculating 
the marginals over momentum or position for each of the 
diagonal terms of the Wigner function, we obtain the 
probability distribution for the conjugate variable. Each 
of the obtained diagonal distributions conserve the initial 
Gaussian shape, both in position and in space, in spite 
of the interaction with the environment, and allow for a 
clear description in terms of the center and width of the 
corresponding Gaussian, which bear a close analogy with 
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the classical Brownian motion of a particle. In particu¬ 
lar, at large times the particle reaches a limit velocity, a 
feature which is particularly important for the descrip¬ 
tion of the SG experiment inside a gas tube. We also 
showed that our results agree with the non decoherence 
expressions in the appropriate limit. 

By adopting realistic parameters for the SG experi¬ 
ment, we plotted the diagonal and off-diagonal elements 
of the Wigner function. We observe that, for low values 
of the decoherence parameter 7 , the initial state sepa¬ 
rates into two components, corresponding to the up and 
down spin, as expected. The off-diagonal terms provide 
us information about decoherence, as a consequence of 
the interaction with the environment, that manifests in 
the damping of these terms as time evolves. In order to 
quantify the effect of decoherence, we evaluate the norm 
of the off-diagonal element over the whole phase space. 
The characteristic decoherence time scale td is defined as 
the instant when the initial norm is reduced by a factor 
e. By obtaining td for different values of 7 , we find a re¬ 
lation which is well described by a power law td oc 7 -1 / 5 . 
Our result differs from previous results, where a depen¬ 
dence td oc 7 1 / 3 was claimed instead. We argue that such 
dependence is counter intuitive, as it would imply that 
a larger value of 7 makes decoherence effects to appear 
later. In contrast, we obtain that a stronger value of 7 
implies a shorter decoherence time, which seems more 
reasonable. 

Finally, we discuss a caveat of the model, that shows 
up after estimating the order of magnitude of 7 in a re¬ 


alistic scenario. In fact, we obtain 7 ~ 10 10 , a value that 
would lead to a large friction, and would ruin the ob¬ 
served phenomenology of the experiment. This estimate 
has itself a loophole, as the calculated mean free path 
of the incoming particles on the tube turns out to be of 
the order of the tube length, thus implying that collisions 
with the air molecules are rare events. Altogether, these 
reasonings rise some doubts on the applicability of the 
Caldeira-Leggett model for this problem. In our opinion, 
this question deserves further research, given the impor¬ 
tance of the experiment as one of the cornerstones in 
quantum physics. 
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V. AUXILIARY FORMULAE FOR THE 

DIAGONAL AND OFF-DIAGONAL WIGNER 
FUNCTION 

In this Appendix, we give the explicit expressions for 
auxiliary functions that define both the diagonal and 
off-diagonal terms of the matrix WF. The functions 
F(z,p,r) and G(r) that enter in the diagonal part Eq. 
(15) are defined as 


F(z,p,r) = 2'yDa 2 ^y 2 m 2 (l — e 2r ) z' 2 — 2'jmp' (l — e T ) 2 z' + p 12 (2 — e r )~+2T + l 

+ 7 4 m 2 p ,2 a 4 + 7 2 h 2 (p' (l — e~ T ) — 7 me~ T z') 


(33) 


G(t) = 8D 2 <j 2 (l — e T ) ((e r + l) r — 2 (l — e T )) + 2^D [ 7 2 to 2 ct 4 (l — e 2r ) 
+ (e" 2 T (2r + 3) - 4e -r + l) h 2 ] + 7 4 m 2 a 2 e~ 2T h 2 , 


(34) 


where z' = z T z c , p' = p TPc- The magnitudes z c and 
p c are given in Eqs. (16) and (17), respectively. 


For the off diagonal term defined in Eq. (18) we intro¬ 
duced the following definitions: 


Ci{t) 


rf A 2 (a 2 (—2 D (2r 3 — 6 r 2 + 6 r + 3) + QDe 2t + 24 De t t — 3^ 3 m 2 a 2 T 2 ) — 3y (1 — r — e T ) 2 h 2 ^j 

3q 5 m 2 <j 2 h 2 


(35) 


C 2 (r) 


ryA ( 7 h 2 (e r — e T r — e 2r ) — 2Da 2 (l — 2e T r — e 2r )) 


7 3 fi 2 TOcr 2 


(36) 
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C 3 (t) = 


a 2 (2 De 2t — 8 De T — D(At — 6 ) — 7 3 ?n 2 tr 2 ) — 7 h 2 (1 — e r ) 2 

47 3 to 2 ct 2 


(37) 


??A ^7 (1 — e T ) (1 — t — e T )h 2 — cr 2 ^2Z) (r + e r — l) 2 + 7 3 m 2 a 2 T^j ^ 
^ 7 4 m 2 a 2 h 


C 5 (t) = 


{e~ T - 1) {^e~ T H 2 + 2Da 2 (1 - e~ r )) 


C 6 (r) = 


27 2 ma 2 h 


7 e 2 r fi 2 + 2D<7 2 (l — e 2t ) 

47 CT 2 fi . 2 


(38) 

(39) 

(40) 
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